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Multiscale  Modeling:  From  Atoms  to  Devices 


2.1  INTRODUCTION 

Coarse-grain  models  provide  an  efficient  means  to  simulate  and  investigate  systems 
in  which  the  desired  behavior,  property,  or  response  is  inherently  at  the  mesoscale — 
those  that  are  both  inaccessible  to  full  atomistic  representations  and  inapplicable 
to  continuum  theory.  Granted,  a  developed  coarse-grain  model  can  only  reflect 
the  behavior  included  in  their  governing  potentials  and  associated  parameters,  and 
consequently,  the  source  of  such  parameters  typically  determines  the  accuracy  and 
utility  of  the  coarse-grain  model.  It  is  our  contention  that  a  complete  theoretical 
foundation  for  any  system  requires  synergistic  multiscale  transitions  from  atomic  to 
mesoscale  to  macroscale  descriptions.  Hierarchical  “handshaking”  at  each  scale  is 
crucial  to  predict  structure -property  relationships,  to  provide  fundamental  mecha¬ 
nistic  understanding,  and  to  enable  predictive  modeling  and  material  optimization  to 
guide  synthetic  design  efforts.  Indeed,  a  finer-trains-coarser  approach  is  not  limited 
to  bridge  atomistic  to  mesoscopic  scales  (which  is  the  focus  of  the  current  discus¬ 
sion),  but  can  also  refer  to  hierarchical  parameterization  transcending  any  scale,  such 
as  mesoscopic  to  continuum  levels.  Such  a  multiscale  modeling  paradigm  establishes 
a  fundamental  link  between  atomistic  behavior  and  the  coarse-grain  representation, 
providing  a  consistent  theoretical  approach  to  develop  coarse-grain  models  for  sys¬ 
tems  of  various  scales,  constituent  materials,  and  intended  applications. 


2.1.1  Motivation:  Hierarchical  Systems  and  Empirical  Links 

Many  biological  tissues  are  composed  of  hierarchical  structures,  which  provide  excep¬ 
tional  mechanical,  optical,  or  chemical  properties  due  to  specific  functional  adapta¬ 
tion  and  optimization  at  all  levels  of  hierarchy.  Nature  has  shown  that  a  material’s 
structure — and  not  its  composition  alone — must  be  considered  in  the  design  of  new 
material  systems  for  use  in  high-performance  applications.  Fundamental  structural 
arrangements  and  the  mechanistic  properties  are  inherently  linked.  Indeed,  nature’s 
integration  of  robustness,  adaptability,  and  multifunctionality  requires  the  merging 
of  structure  and  material  across  a  broad  range  of  length  scales,  from  nano  to  macro, 
and  is  apparent  in  biological  materials  such  as  bone,  wood,  and  protein-based  materi¬ 
als  [1-3].  The  analysis  of  such  hierarchical  materials  is  an  emerging  field  that  uses 
the  relationships  between  multiscale  structures,  processes,  and  properties  to  probe 
deformation  and  failure  phenomena  at  the  molecular  and  microscopic  levels  [4]. 

For  the  current  discussion,  the  term  hierarchical  is  used  loosely  to  indicate  a 
material  system  with  at  least  a  single  distinct  differentiation  between  constituent 
material  components  and  global  system  structure.  For  example,  in  Chapter  3  we  dis¬ 
cuss  both  coarse-grain  modeling  of  carbon  nanotube  arrays  and  collagen  fibrils.  For 
the  nanotube  arrays,  the  components  are  defined  by  individual  carbon  nanotubes, 
and  thus  the  array  is  considered  a  hierarchical  structure.  For  the  collagen  fibrils, 
the  components  are  defined  as  tropocollagen  molecules,  while  the  system  of  inter¬ 
est  is  the  entire  fibril.  It  is  noted  that  tropocollagen  fibrils  are  themselves  composed 
of  a  hierarchical  arrangement  of  polypeptide  chains,  which  are  also  composed  of 
constituent  amino  acids.  Thus,  the  defined  components  need  not  be  the  fundamen¬ 
tal  building  blocks  of  the  system.  In  contrast,  we  define  the  components  of  alpha- 
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helical  proteins  (see  case  study  in  Chapter  3)  as  a  single-protein  convolution,  while 
the  system  is  characterized  by  the  entire  protein,  recognizing  the  hierarchical  effects 
of  single  molecular  conformations.  The  coarse-graining  procedures  discussed  here 
focus  on  constituent  materials  that  can  be  modeled  by  full  atomistic  techniques  using 
classical  molecular  dynamics,  while  the  system  structure  (and  relevant  behaviors) 
requires  a  fully  informed  coarse-grain  representation.  The  atomistic  model  provides 
the  underlying  physics  to  the  coarse-grain  potentials,  allowing  the  investigation  of 
structure-property  functions  at  the  required  length  scale. 

As  an  analogy,  such  multiscale  model  transitions  can  be  compared  to  the  parame¬ 
terization  of  atomistic  force  fields  by  detailed  quantum  mechanical  methods  (Figure 
2.1).  It  is  possible  to  derive  the  force  constant  parameters  used  in  full  atomistic  mod¬ 
els  from  first  principle  quantum  mechanical  calculations  [5].  More  sophisticated 
techniques,  such  as  density  functional  theory  (DFT),  are  simplified  in  the  parameters 
and  functions  of  atomistic  potentials.  The  atomic  structure  (bosons,  fermions,  etc.) 
can  be  considered  as  the  material  components,  while  bonded  atoms  (such  as  carbon- 
carbon  bonds)  can  be  thought  of  as  the  structural  system.  From  this  perspective,  atom¬ 
istic  force  fields  can  be  considered  as  a  thorough  and  systematic  coarse-graining  of 
quantum  mechanics.  The  direct  application  of  quantum  mechanical  results  to  large 
macromolecular  systems  (such  as  proteins)  is  not  yet  feasible  because  of  the  large 
number  of  atoms  involved.  As  a  logical  extension,  we  analyze  the  results  of  full 
atomistic  simulations  and  implement  them  in  the  fitting  of  mesoscopic  coarse-grain 
potentials.  Each  transition  involves  an  increase  in  accessible  system  size,  length  scale, 
and  time  scale,  as  well  as  a  simplification  of  the  governing  theoretical  models.  As  a 
consequence,  details  of  interactions  are  lost,  but  the  desired  fundamental  behavior 


Time  and  length  scales 

FIGURE  2.1  Rationale  for  coarse-graining  approach  through  analogous  comparison 
between  the  parameterization  of  atomistic  force  fields  by  detailed  quantum  mechanical  meth¬ 
ods  to  the  development  of  coarse-grain  potentials  by  the  results  of  full  atomistic  simulations, 
(a)  Electron  charge  isosurface  plot  of  (5,5)  carbon  nanotube  via  density  functional  theory 
(DFT).  (b)  Segment  of  full  atomistic  model  of  identical  carbon  nanotube  using  carbon- 
carbon  interactions,  (c)  Top  view  of  adhered  bundle  of  coarse-grain  carbon  nanotubes  using 
potentials  developed  from  atomistic  results.  Each  transition  involves  an  increase  in  accessible 
system  size,  length  scale,  and  time  scale,  as  well  as  a  simplification  of  the  governing  theoreti¬ 
cal  models. 
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is  maintained.  It  is  noted  that  this  “desired  fundamental  behavior”  is  dependent  on 
the  system  application,  and  can  include  mechanistic  behaviors,  thermal  or  electrical 
properties,  molecular  interactions,  equilibrium  configurations,  etc.  To  illustrate,  an 
atomistic  representation  of  carbon  bonding  can  include  terms  for  atom  separation, 
angle,  and  bond  order  (such  as  CHARMM  [6]-  or  AMBER  [6,7]-type  force  fields), 
but  lacks  any  description  of  electron  structure,  band  gaps,  or  transition  states  found  in 
quantum  mechanical  approaches  such  as  DFT.  However,  the  atomistic  model  main¬ 
tains  the  accurate  behavior  of  the  individual  carbon  atoms  and  neglecting  the  effects 
of  electrons  is  deemed  a  necessary  simplification.  Likewise,  the  coarse-graining 
of  a  carbon  nanotube  integrates  the  effects  of  multiple  carbon  bonds  into  a  single 
potential.  The  exact  distribution  of  carbon  interactions  is  lost,  but  the  behavior  of  the 
molecular  structure  is  maintained. 

The  development  of  coarse-grain  models  allows  the  simulation  of  events  on  physi¬ 
cal  time  and  size  scales,  leading  to  a  range  of  possible  advances  in  nanoscale  design 
and  molecular  engineering  in  a  completely  integrated  bottom-up  approach.  A  coarse- 
graining  approach  is  intended  to  develop  tools  to  investigate  material  properties  and 
underlying  mechanical  behavior  typically  required  for  material  design  that  system¬ 
atically  integrates  characteristic  chemical  responses.  However,  it  is  emphasized  that 
the  intention  is  not  to  circumvent  the  need  for  full  atomistic  simulations — in  con¬ 
trast,  coarse-graining  requires  accurate  full  atomistic  representations  to  acquire  the 
necessary  potential  parameters.  The  hner-trains-coarser  procedure  described  here 
is  an  attempt  to  reconcile  first  principles  derivations  with  hierarchical  multiscale 
techniques  by  using  atomistic  theory  with  molecular  dynamics  simulations  in  lieu  of 
empirical  observations  in  a  unified  and  systematic  approach. 


2.1.2  Diversity  of  Systems  and  Applications: 

Need  for  a  System-Dependent  Approach 

Atomistic  force  fields,  which  must  encompass  atom-atom  bonds  and  interactions, 
can  be  developed  in  a  general  formulation  due  to  the  common  molecular  compo¬ 
nents.  For  example,  the  behavior  of  carbon-carbon  bonds  in  a  carbon  nanotube  is 
similar  to  the  alpha-carbon  backbone  of  a  protein  sequence  in  terms  of  strength  and 
bond  length.  A  thoroughly  developed  atomistic  force  field  is  capable  of  representing 
a  vast  assortment  of  molecular  systems.  Coarse-grained  potentials,  conversely,  are 
usually  developed  to  represent  a  particular  system,  and  consequently  feature  unique 
idiosyncrasies  in  construction  and  behavior.  Accordingly,  differences  in  system  com¬ 
plexity  and  goals  of  modeling  lead  to  difficulties  in  developing  a  universal  method 
for  coarse-graining. 

Attempts  to  avoid  this  system-dependent  approach  (i.e.,  generalized  coarse- 
graining  frameworks)  can  potentially  result  in  a  complex  coarse-grain  description  to 
account  for  the  multitude  of  molecular  interactions  a  general  description  must  encom¬ 
pass.  The  introduction  of  many  potentials  and  parameters  essentially  mimic  the 
function  and  form  of  full  atomistic  force  fields,  albeit  at  a  coarse-grain  scale.  Other 
generalized  frameworks  attempt  to  simplify  the  description  of  molecular  interactions 
as  much  as  possible.  For  example  the  MARTINI  force  field  [8]  has  been  successful 
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in  the  intent  of  efficiently  modeling  larger  systems,  and  warrants  further  discussion 
(see  Section  2.3).  To  capture  fundamental  interactions,  general  coarse-grain  poten¬ 
tials  developed  for  application  to  multiple  systems  are  limited  to  few  atom  mappings 
(i.e.,  two-,  four-,  and  six-bead  models)  and  are  beneficial  and  appropriate  for  systems 
where  the  interactions  are  still  at  the  atomistic  scale  (i.e.,  proteins,  polymers,  etc.). 

It  is  apparent  that  different  material  systems  can  be  characterized  by  mechanical 
properties  at  the  atomistic  to  microscale,  requiring  a  more  general  coarse-graining 
approach  to  model  the  intended  structure -property  behavior.  It  is  our  proposition 
that  a  hierarchical  multiscale  system  requires  a  system-dependent  approach  to  coarse- 
graining.  Accurate  system  representation  is  maintained  at  a  cost  of  coarse-grain 
potential  versatility.  By  applying  a  system-dependent  finer-trains-coarser  approach, 
a  variety  of  systems  can  be  coarse-grained  and  modeled  for  different  purposes,  tran¬ 
scending  different  scales  and  functionalities.  Such  derived  coarse-grain  potentials 
are  not  meant  to  be  universally  applicable,  but  provide  a  means  to  investigate  a  spe¬ 
cific  system  under  specific  conditions. 


2.1.3  Alternative  Coarse-Graining  Advantages: 

Pragmatic  System  Simplification 

A  commonly  stated  motivation  and  presumed  primary  benefit  for  the  development 
of  coarse-grain  potentials  and  models  is  to  allow  the  simulation  of  larger  systems 
at  longer  time  scales.  Indeed,  the  reduction  of  system  degrees-of-freedom  and  the 
smoother  potentials  implemented  allow  larger  time-step  increments  (and  thus  time 
scales)  and  each  element  is  typically  an  order  of  magnitude  or  more  larger  in  length 
scale.  Additionally,  cheaper  potential  calculations  (in  terms  of  computational  effi¬ 
ciency)  can  be  exploited  to  either  increase  the  number  of  coarse-grain  elements  (thus 
representing  even  larger  systems)  or  simulate  a  relatively  small  system  over  more 
integration  steps  (further  extending  accessible  time  scales).  Such  benefits  are  inher¬ 
ent  to  any  coarse-grain  representation,  and  can  serve  as  a  de  facto  definition  of  the 
coarse-graining  approach. 

Nevertheless,  a  pragmatic  approach  of  system  simplification  is  found  in  many 
engineering  disciplines.  Complex  electronic  components  are  designed  based  on  sim¬ 
plified  models  of  circuits,  with  element  behavior  defined  by  such  general  proper¬ 
ties  as  current,  voltage,  and  resistance.  Robust  building  structures  are  analyzed  via 
notions  of  beam  deflections  and  beam-column  joint  rotations,  among  other  simpli¬ 
fying  assumptions.  In  both  cases,  more  detailed  system  representations  are  known 
and  can  be  implemented  (e.g.,  implementation  of  temperature  and  material  effects 
in  a  transistor,  or  a  detailed  frame  analysis  including  stress  concentrations  in  bolts 
or  welds).  It  is  apparent  that  such  additions  result  in  a  more  accurate  representa¬ 
tion  of  the  modeled  system,  but  also  serve  to  increase  the  computational  expense 
of  analysis  as  well  as  introduce  a  more  sophisticated  theoretical  framework  (which 
subsequently  requires  a  more  detailed  set  of  material  and  model  parameters).  Rarely, 
however,  is  the  use  of  simplified  and  computationally  efficient  models  justified  by 
inaccessible  time  and  length  scales  of  the  more  detailed  description.  Such  models 
are  applied  for  analysis  in  lieu  of  a  more  detailed  description  because  they  provide 
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an  accurate  representation  of  the  system-level  behavior  and  response  with  confidence 
in  the  properties  of  the  model  components. 

For  example,  a  common  application  of  system-level  analysis  is  used  in  the  study  of 
structural  frames.  Required  loads  and  resistances  are  computed  at  the  system-level, 
and  the  structural  components  are  assumed  to  behave  according  to  a  generalized  theo¬ 
retical  framework  (in  terms  of  beam  deflections,  rotations,  etc.).  Detailed  structural 
elements,  such  as  steel  joist  girders,  are  presumed  to  behave  according  to  their  repre¬ 
sentative  models  and  do  not  warrant  a  full  analysis  at  the  structural  level.  Of  course,  the 
joist  girders  themselves  are  rigorously  tested  at  the  component  level  to  determine  the 
necessary  parameters  (ultimate  loads,  yield  stresses,  etc.).  We  apply  the  same  approach 
to  molecular  coarse-graining.  Full  atomistic  simulations  are  used  to  probe  and  acquire 
the  behavior  of  molecular  components  and  parameterize  coarse-grain  potentials.  These 
potentials,  in  turn,  are  implemented  in  coarse-grain  representations  that  can  facilitate 
the  analysis  of  large-scale  system-level  phenomenon.  There  is  an  inherent  assumption 
that  the  coarse-grain  theoretical  framework  provides  an  accurate  representation  of  the 
atomistic  behavior,  just  as  application  of  elastic  beam  theory  assumes  the  reliability  of 
a  steel  joist  girder.  See  Figure  2.2  for  a  schematic  comparison. 

It  is  stressed  that  even  if  full  atomistic  representations  are  computationally  possi¬ 
ble,  a  coarse-grain  description  can  still  be  suitable  for  systematic  analysis  of  variable 
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Detailed  model  of  constituent 
material(s)  and  structures(s) 


Steel  joist  girder  with  known 
member  properties  (Yield  stress, 
Young’s  modulus,  etc.)  and 
geometry. 


Full  atomistic  model  of  amyloid  with 
known  atomistic  interactions  and 
chemical  structure. 


Joist  girder  represented  by  a  beam 
element  with  known  behavior  and  load 
response  from  analysis  of  joist  girder; 
assumed  linear  elastic  beam  theory 
valid. 


Amyloid  represented  by  mesoscale 
model  elements  with  coarse-grain 
potentials  developed  based  on  known 
behavior  and  response  of  full  atomistic 
simulations. 


Simplified  system-level 
analysis 


Beam  elements  used  in  analysis  of 
structural  frame. 


Coarse-grain  elements  used  in 
analysis  of  amyloid  plaque 
(entanglement  of  many  fibers),  a 
protein  deposit  associated  with 
Alzheimer’s  disease. 


FIGURE  2.2  (See  color  insert  following  page  146.)  Analogous  comparison  of  system 
simplification  between  structural  analysis  and  materiomics.  Detailed  model  constituents 
and  their  structural  arrangement  are  analyzed  to  parameterize  a  coarse-grain  representa¬ 
tive  model,  maintaining  relevant  behaviors  and  implemented  for  simplified  analysis.  Here, 
we  see  the  transition  of  steel  joist  girders  to  beam  elements  to  structural  frame  analysis, 
paralleled  by  a  full  atomistic  model  of  an  amyloid,  a  coarse-grain  mesoscopic  representa¬ 
tion,  and  the  system-level  amyloid  plaque  (aggregation  of  thousands  of  amyloids).  (Image  of 
amyloid  plaque  reprinted  with  permission  from  Macmillan  Publishers  Ltd.,  Pilcher,  H.R., 
Nature  News,  2003.  Copyright  2003.) 
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system  configurations,  which  requires  a  large  number  of  simulations.  We  offer  an 
alternative  motivation  for  a  coarse-graining  approach  complementary  to  extension 
of  accessible  time  and  length  scales,  where  the  catalyst  for  coarse-grain  potential 
development  is  not  the  extension  of  traditional  molecular  dynamics,  but  to  provide 
an  accurate  and  reliable  method  for  system-level  analysis  and  probe  the  mechani¬ 
cal  response  and  structure-property  relation  for  hierarchical  systems.  Furthermore, 
mesoscopic  models  provide  simple  and  efficient  modeling  techniques  for  experimen¬ 
talists,  allowing  a  more  direct  comparison  between  simulation  and  a  vast  variety  of 
experimental  techniques,  without  requiring  specialized  molecular  dynamics  tools 
such  as  specialized  computer  clusters  with  complex  software. 

2.1.4  When  to  Coarse-Grain:  Appropriate  Systems  and  Considerations 

It  is  emphasized  that  not  all  systems  will  benefit  from  a  coarse-grain  representa¬ 
tion  and  prudent  consideration  must  be  given  regarding  the  system  characterization 
and  intent  of  the  simulations.  For  some  systems,  complex  behaviors  may  require 
full  atomistic  representation,  or  material  inhomogeneities  may  not  be  able  to  be 
described  by  coarse-grain  elements.  Such  systems  may  benefit  and  indeed  require 
full  atomistic  representations.  Typical  motivations  for  a  coarse-graining  approach 
include: 

1.  Inaccessible  time  scale  for  phenomenon  or  behavior  via  full  atomistic 
representation 

2.  Inaccessible  length  scale  for  phenomenon  or  behavior  via  full  atomistic 
representation 

3.  Focus  on  global  system  properties  and/or  mechanical  behavior  rather  than 
on  molecular  structure  and/or  chemical  interactions 

4.  Desire  for  a  direct  simplified  analysis  of  simulation  results  and  system 
behavior 

In  addition  to  these  motivating  factors,  some  systems  are  more  conducive  to 
coarse-graining  due  to  chemical  composure  and/or  structural  arrangement.  Coarse- 
grain  models  are  easily  adapted  for  atomistically  homogenous  systems,  consisting  of 
repetitive  structures  such  as  carbon  nanotubes  (with  a  uniform  cylindrical  nanostruc¬ 
ture).  By  extension,  atomistically  heterogeneous  materials,  such  as  protein-based 
materials,  are  applicable  to  coarse-graining  if  they  are  mesoscopically  homogeneous, 
where  the  local  effects  of  distinct  amino  acids  (or  other  molecular  inhomogeneities) 
produce  similar  global  behaviors  and  are  deemed  negligible.  Other  system  properties 
to  consider  include  the  discretization  or  material  units  and/or  mechanical  behavior, 
with  a  logical  correlation  to  coarse-grain  elements,  and  any  hierarchical  structures  in 
which  the  intended  structure-property  relation  is  to  be  investigated. 

2.2  EXAMPLES  OF  COARSE-GRAINING  METHODS 

In  the  past  decade,  various  simple  models  have  been  used  to  describe  the  large- 
scale  motions  of  complex  molecular  structures  where  more  detailed  classical 
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phenomenological  potentials  [6,7]  involving  all  atoms  cannot  be  used  because  of  the 
restrictions  on  the  amount  of  time  that  can  be  covered  in  computer  simulations.  The 
reader  is  referred  to  recent  reviews  for  a  more  thorough  discussion  of  techniques  and 
applications  [9-12]. 

Single -bead  models  are  the  most  direct  approach  taken  for  studying  macromol¬ 
ecules.  The  term  single  bead  derives  from  the  idea  of  using  single  beads,  that  is, 
point  masses,  for  describing  a  functional  group,  such  as  single  amino  acids,  in  a  mac- 
romolecular  structure.  The  elastic  network  model  (ENM)  [13],  Gaussian  network 
model  [14],  and  GO  model  [15]  are  well-known  examples  that  are  based  on  such  bead 
model  approximations.  These  models  treat  each  amino  acid  as  a  single  bead  located 
at  the  Ca  position,  with  mass  equal  to  the  mass  of  the  amino  acid.  The  beads  are 
connected  via  harmonic  bonding  potentials,  which  represent  the  covalently  bonded 
protein  backbone.  Elastic  network  models  have  been  used  to  study  the  properties 
of  coarse-grained  models  of  proteins  and  larger  biomolecular  complexes,  focusing 
on  the  structural  fluctuations  about  a  prescribed  equilibrium  configuration,  such  as 
normal  vibrational  modes,  and  are  discussed  further  in  Section  2.2.1. 

In  GO-like  models,  an  additional  Lennard-Jones-based  term  is  included  in  the 
potential  to  describe  short-range  nonbonded  interactions  between  atoms  within  a 
finite  cutoff  separation.  Despite  their  simplicity,  these  models  have  been  extremely 
successful  in  explaining  thermal  fluctuations  of  proteins  [11]  and  have  also  been 
implemented  to  model  the  unfolding  problem  to  elucidate  atomic -level  details  of 
deformation  and  rupture  that  complement  experimental  results  [16-18].  A  more 
recent  direction  is  coupling  of  ENM  models  with  a  finite  element-type  framework 
for  mechanistic  studies  of  protein  structures  and  assemblies  [19]. 

Using  more  than  one  bead  per  amino  acid  provides  a  more  sophisticated  descrip¬ 
tion  of  protein  molecules.  In  the  simplest  case,  the  addition  of  another  bead  can  be 
used  to  describe  specific  side-chain  interactions  in  proteins  [20].  Higher-level  models 
(e.g.,  four-  to  six-bead  descriptions)  capture  more  details  by  explicit  or  united  atom 
description  for  backbone  carbon  atoms,  side  chains,  and  carboxyl  and  amino  groups 
of  amino  acids  [21,22],  Coarser-level  multiscale  modeling  methods  similar  to  those 
presented  here  have  been  reported  more  recently,  applied  to  model  biomolecular 
systems  at  larger  time  and  length  scales.  These  models  typically  employ  super-atom 
descriptions  that  treat  clusters  of  amino  acids  as  beads.  In  such  models,  the  elasticity 
of  a  polypeptide  chain  is  captured  by  simple  harmonic  or  anharmonic  (nonlinear) 
bond  and  angle  terms.  These  methods  are  computationally  quite  efficient  and  cap¬ 
ture  shape-dependent  mechanical  phenomena  in  large  biomolecular  structures  [23], 
and  can  also  be  applied  to  collagen  fibrils  in  connective  tissue  [24]  as  well  as  min¬ 
eralized  composites  such  as  nascent  bone  [25].  The  development  of  the  coarse-grain 
model  for  collagen  will  be  discussed  further  as  a  case  study  in  Chapter  3. 

We  proceed  to  provide  a  brief  discussion  of  some  of  the  aforementioned  approaches 
implemented  for  coarse-grain  potentials  as  an  overview  of  other  methods.  The  fol¬ 
lowing  discussion  is  not  intended  to  include  all  the  intricacies  and  details  of  the 
development  and  application  of  each  method,  but  merely  to  provide  examples  of 
various  types  of  coarse-graining  procedures  and  the  progression  of  adding  more 
complexity  to  the  coarse-grain  representation. 
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2.2.1  Elastic  Network  Models 

In  the  simplest  form,  a  coarse-grain  model  can  be  defined  by  a  single  potential  for 
all  beads.  For  example,  the  aforementioned  ENMs  can  be  thought  of  as  a  single  pair 
potential  between  neighbor  atoms,  such  that: 

^CG  =  ^ENM  =  ^^0ENM(r)>  (2.1) 

where 

0ENM(r)  =  ~2^r  {r  ~  rfl)  (2-2) 

Here,  each  pair  of  atoms  is  assigned  a  harmonic  spring  bond  potential  with  stiffness, 
Kr,  about  an  initial  equilibrium  spacing,  r0  (see  Figure  2.3  for  an  example). 

The  elastic  constant,  Kr,  can  be  defined  as  a  constant  [13]  or,  as  another  example, 
as  a  function  of  r0\ 


K,{r,)  =  K  !- 
\r0  J 


(2.3) 


(b)  Elastic  network  model 


FIGURE  2.3  (See  color  insert  following  page  146.)  (a)  Full  atomistic  representation  and 
(b)  elastic  spring  network  representation  of  an  amyloid  fibril.  Here,  springs  are  connected 
to  all  neighbor  atoms  within  a  cutoff  of  10  A,  and  spring  stiffnesses  are  assigned  accord¬ 
ing  to  an  exponentially  decaying  function  (Equation  2.4).  Model  was  implemented  to  deter¬ 
mine  normal  vibration  modes  and  structure  stiffness.  (Model  images  courtesy  of  Dr.  Z.  Xu, 
Massachusetts  Institute  of  Technology.) 
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Kr(r0)  =  Kex  p 


r  \2 

>b 


(2.4) 


Equation  2.3  scales  the  mean  value  of  the  elastic  constant,  K,  with  the  initial  bond 
length,  r0,  such  that  Kr  ( r  j  =  K.  If  the  ENM  is  representative  of  a  solid  or  continuous 
media,  this  method  treats  all  elastic  elements  as  if  they  had  the  same  cross-sectional 
area.  A,  and  a  constant  Young’s  modulus,  E,  such  that  the  quantity  (EA/r0)  is  the 
same  for  all  elements  in  the  network  [26].  The  second  method  subjects  the  elastic 
constant  to  an  exponentially  decaying  function  with  respect  to  a  mean  bond  distance 
(r).  Such  a  formulation  of  Kr  can  represent  weak  interactions  of  atoms  at  a  distance 
(such  as  van  der  Waals  interactions)  and  provide  a  more  complex  description  of  inter¬ 
actions  [27,28]. 

Such  ENMs  fall  into  two  broad  classes,  those  with  Hookean  springs  that  describe 
the  rigidity  of  the  macromolecule  [13],  and  those  describing  the  connectivity  of  the 
macromolecule  [29].  For  illustrative  purposes,  we  limit  our  discussion  to  scalar 
spring  constants  and  forces  along  the  vector  between  two  atoms,  but  also  note  that 
more  complex  model  representations  introducing  directionality  and  anisotropy  can 
be  formulated  [30].  As  such,  even  a  single  potential  description  can  become  increas¬ 
ingly  sophisticated  as  applications  attempt  to  probe  complex  molecular  deformations 
such  as  protein  residue  fluctuations  [31]  and  equilibrium  state  transitions  [32,33]. 


2.2.2  Two  Potential  Freely  Jointed  Chain  Polymer  Models 

Although  useful  in  the  application  of  normal  mode  analysis  of  single  protein  macro¬ 
molecules,  ENM  techniques  lack  the  description  necessary  for  intermolecular  inter¬ 
actions  and  large  deformation  from  equilibrium  conditions.  The  subsequent  step 
in  coarse-grain  model  development  is  the  combination  of  two  simple  potentials  to 
represent  the  intermolecular  and  intramolecular  interactions  of  individual  macro¬ 
molecules.  Polymer  systems  are  frequently  represented  by  two  such  coarse-grain 
potentials,  which  encompass  intramolecular  (bonded)  interactions  and  intermolecu¬ 
lar  (nonbonded)  interactions,  respectively,  in  a  freely  jointed  chain  (FJC)  representa¬ 
tion  (no  angular  constraints)  or: 

^CG  =  ^bonded  +  ^nonbonded  =  0FENE  (r)  +  0WCA  (r)  (2-5) 

For  bonded  beads,  a  finitely  extensible  nonlinear  elastic  (FENE)  potential  [34-36] 
is  implemented  to  maintain  distance  between  connected  beads  and  prevent  polymer 
chains  from  crossing  each  other: 


0fene  (r)  —  „  kr0  In 


( 

2~ 

1- 

r 

Jo) 

for  r  <  r0> 


(2.6) 


4>feneW  =  00  for  r  >  r0 
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The  Weeks-Chandler-Anderson  (WCA)  potential  [37,38]  is  the  Lennard-Jones 
12:6  potential  truncated  at  the  position  of  the  minimum  and  shifted  to  eliminate 
discontinuity: 


0wca  (r)  =  4e 


f  ^ 

G 

\r  J 


f  \ 

G 

\V  J 


for  —  <  \f2, 
a 


(2.7) 


0wcaW  =  O  for  —  >^2 

a 

The  WCA  potential  results  in  a  purely  repulsive  potential.  The  combination  of  the 
FENE  attractive  potential  with  the  WCA  repulsive  interaction  creates  a  potential 
well  for  the  flexible  bonds  that  can  maintain  the  topology  of  the  molecule  [39].  Such 
models  can  successfully  represent  stretching,  orientation,  and  deformation  of  poly¬ 
mer  chains  and  simple  biomolecules. 

Due  to  the  inherent  flexibility  of  the  freely  jointed  chains,  this  coarse-graining 
approach  is  particularly  suited  for  systems  defined  by  long-chain  polymers  with 
relatively  short  persistence  lengths,  or  systems  that  are  entropically  driven.  As  such, 
example  applications  include  the  investigation  of  viscoelastic  behavior  of  polymer 
melts  [40],  stretching  of  polymers  in  flow  [41],  and  other  rheological  properties 
[42,43].  More  complex  formulations  use  the  combination  of  FENE  and  WCA  poten¬ 
tials  to  define  the  coarse-grain  model,  and  add  additional  potentials  for  a  more  robust 
description  of  the  system,  such  as  the  introduction  of  angular  terms  (for  stiffness) 
and  Coulombic  terms  (for  electrostatic  interactions)  to  model  complex  macromol¬ 
ecules  such  as  DNA  [44]. 

2.2.3  Generalization  of  Interactions:  The  MARTINI  Force  Field 

The  lack  of  intermolecular  interaction  characterization  in  the  previous  coarse- 
graining  approaches  required  a  more  complex  formulation — one  that  maintained 
a  definitive  structure  of  the  macromolecule,  as  well  as  integrated  the  interactions 
between  macromolecules.  The  development  of  the  MARTINI  force  field  attempted 
to  provide  a  general  coarse-grain  model  that  could  be  efficiently  adapted  for  a  mul¬ 
titude  of  biological  systems  by  taking  advantage  of  the  fact  that  the  majority  of  bio¬ 
logical  molecules  (such  as  protein  structures  and  lipids)  are  composed  of  the  same 
categories  of  functional  groups  at  the  atomistic  level  [8,45,46].  Such  approaches 
were  previously  implemented  with  various  degrees  of  complexity  for  application  to 
specific  systems,  integrating  multiple  beads  per  functional  group,  or  more  complex 
parameter  formulations  (see  Marrink  et  al.  [45]  for  a  discussion  and  Shelley  et  al. 
[47,48]  as  examples).  Essentially,  by  coarse-graining  the  functional  groups,  mole¬ 
cules  with  different  architectures  can  be  easily  built  and  simulated.  The  aim  of  the 
MARTINI  force  field  was  to  retain  the  chemical  nature  of  the  molecular  components 
while  defining  as  few  bead  types  as  possible.  The  applied  philosophy  was  to  avoid 
a  focus  on  the  reproduction  of  structural  details  for  a  particular  system,  but  rather 
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aim  for  a  broad  range  of  applications  without  the  need  to  reparameterize  the  model 
each  time.  As  a  result,  there  is  a  slight  trade-off  between  accuracy  and  applicability, 
with  many  different  biological  systems  able  to  be  modeled.  Only  four  interaction 
types  are  considered — polar,  nonpolar,  apolar,  and  charged — and  subtypes  of  these 
interactions  are  used  to  define  hydrogen-bonding  characteristics  (see  Figure  2.4  for 
coarse-grain  representations  of  amino  acids). 

The  MARTINI  force  field  can  be  represented  by 

^MARTINI  =  ^bond  +  ^angle  +  ^nonbonded  ~  0T  (r)  +  Qe  W  +  ‘Pu  (r)’  (2-8) 

where  harmonic  potentials  are  implemented  for  the  bonded  and  angle  potentials,  or: 

<h(r)=^kT(r-r0f  (2.9) 


Mr) 


kg  (cos£>-  cos 


(2.10) 


The  harmonic  potentials  are  parameterized  by  constants  ( kT ,  r0,  ke,  and  0{])  common 
to  all  bead-types.  Such  an  approach  increases  the  general  applicability  and  versatil¬ 
ity  of  the  force  field.  Nonbonded  interactions  between  interaction  sites  i  and  j  are 
described  by  the  Lennard-Jones  12:6  potential: 
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FIGURE  2.4  (See  color  insert  following  page  146.)  Coarse-grain  representations  of  all 
amino  acid  types  for  MARTINI  force  field  formulation  for  proteins.  Different  colors  repre¬ 
sent  different  particle  types  consisting  of  four  main  types  of  interaction  sites:  polar,  nonpolar, 
apolar,  and  charged.  (Reprinted  with  permission  from  Monticelli,  L.,  et  al.,  J.  Chem.  Theory 
Comput.,  4,  819-834,  2008.  Copyright  2008  American  Chemical  Society.) 
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with  Ojj  representing  the  effective  minimum  distance  between  two  particles,  and 
the  strength  of  their  interaction.  The  uniqueness  of  the  MARTINI  force  field  param¬ 
eterization  is  the  arrangement  of  both  <rI(  and  etj  into  distinct  subsets  of  interaction 
groups. 

The  major  advantage  of  the  MARTINI  force  field  is  the  simplicity  and  versa¬ 
tility  of  the  parameterization — many  proteins,  lipids,  and  molecules  can  be  built 
using  the  same  set  of  coarse-grain  building  blocks.  However,  these  MARTINI-type 
formulations  restrict  the  coarse-graining  to  distinct  functional  groups  to  accurately 
maintain  atomistic  interactions  while  increasing  computational  efficiency.  The  obvi¬ 
ous  disadvantage  of  such  approaches  lies  in  the  restriction  of  the  coarse-graining 
scale.  For  example,  on  average,  each  bead  of  the  MARTINI  model  has  the  volume 
of  four  water  molecules.  As  a  result,  the  time  and  length  scales  of  MARTINI  and 
MARTINI-type  simulations  are  still  limited,  albeit  much  more  efficient  than  full 
atomistic  approaches.  Again,  prudence  is  required  in  selecting  the  level  of  detail  of 
the  coarse-grain  representation. 

For  appropriate  systems  such  as  biological  membranes  consisting  of  lipid  bilay¬ 
ers  (Figure  2.5)  and  lipid-protein  interactions,  such  approaches  provide  a  powerful 
simulation  framework  (see  Venturoli  et  al.  [49]  and  Nielson  et  al.  [50]  for  more  thor¬ 
ough  reviews  and  discussion  of  such  methods). 


FIGURE  2.5  (See  color  insert  following  page  146.)  Simulation  snapshot  of  DPPC/choles- 
terol  bilayer  structure  with  MARTINI  coarse-grain  model  representation.  Cholesterol  mol¬ 
ecules  are  displayed  in  green,  with  red  hydroxyl  groups.  The  DPPC  lipid  tails  are  shown  in 
silver.  Lipid  head  groups  are  displayed  in  purple  and  blue.  Such  a  complex  system  can  only 
be  simulated  via  a  coarse-graining  approach.  (Reprinted  with  permission  from  Marrink,  S.J., 
de  Vries,  A.H.,  and  Mark,  A.E.,  J.  Phys.  Chem.  B,  108,  750-760,  2004.  Copyright  2007 
American  Chemical  Society.) 
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2.2.4  Universal  Framework,  Diverse  Applications 

The  aforementioned  coarse-graining  techniques,  ENM  models,  FJC  representations, 
and  MARTINI-type  formulations  have  been  successful  in  their  own  applications. 
ENMs  reflect  the  atomistic  structure,  FJC  representations  can  accurately  reproduce 
polymer  flow  and  rheological  properties,  while  MARTINI-type  models  accurately 
capture  the  interactions  of  functional  groups  for  applications  in  assembly.  However, 
such  techniques  focus  on  molecular  geometry  and  interactions  rather  than  reflect 
accurate  mechanical  properties  and  response.  We  predict  that  the  investigation  of  the 
mechanical  behavior  of  hierarchical  systems  requires  a  system-dependent  param¬ 
eterization  of  coarse-grain  potentials,  with  a  focus  on  maintaining  both  molecular 
interactions  and  molecular  mechanics.  Thus,  we  introduce  a  universal  framework 
through  a  finer-trains-coarser  multiscale  paradigm,  which  effectively  defines  coarse- 
grain  potentials  via  the  response  of  full  atomistic  simulations,  while  introducing 
relevant  mechanical  properties  at  the  mesoscopic  scale.  The  approach  is  unique  in 
its  ability  to  transcend  multiple  scales,  enabling  the  investigation  of  a  broad  regime 
of  systems  from  protein-based  materials  to  synthetic  composite  structures,  while 
emphasizing  a  theoretical  foundation  on  full  atomistic  investigations  and  asserting 
energetic  equivalence  and  consistent  mechanical  behavior  between  all  levels  of  mod¬ 
eling.  The  result  is  a  set  of  problem-specific  coarse-grain  representations  with  diverse 
applications.  The  discussion  in  Section  2.3  presents  guiding  principles  behind  model 
formulation. 

2.3  MODEL  FORMULATION 

From  the  previous  discussion,  it  is  apparent  that  a  universal,  stepwise  procedure 
for  coarse-graining  a  system  is  limited  by  the  inherent  simplifications  and  intent  of 
the  coarse-grain  representation.  A  model  developed  to  investigate  arrays  of  carbon 
nanotubes  (see  Chapter  3,  Case  Study  I)  will  differ  from  a  model  developed  for  the 
simulation  of  alpha-helix  unfolding  (see  Chapter  3,  Case  Study  II),  even  though  the 
coarse-grain  elements  and  potentials  are  similar.  The  coarse-graining  of  a  particular 
system  can  be  characterized  by  trial-and-error  and  subjective  omissions  or  inclu¬ 
sions  of  pertinent  behaviors,  and  as  such  considered  just  as  much  an  art  as  a  science. 
Careful  consideration  must  be  taken  for  the  intended  application  and  purpose  of  the 
coarse-grain  model.  Indeed,  such  a  judicious  approach  is  the  motivation  for  a  finer- 
trains-coarser  approach.  A  deliberate  focus  on  the  atomistic  behavior  to  characterize 
the  coarse-grain  model  not  only  provides  a  complete  “bottom-up”  theoretical  basis 
for  a  material  system,  but  also  assists  in  the  formulation  of  the  coarse-grain  model 
by  delineating  relevant  behavior  and  interactions.  Nevertheless,  although  we  con¬ 
tend  that  a  step-by-step  coarse-graining  “recipe”  is  problematic,  we  can  outline  the 
“ingredients”  in  a  general,  coarse-graining  framework: 

1.  Potentials  required  to  define  and  characterize  the  mesoscale  system 

2.  Full  atomistic  “test  suite”  applied  to  determine  required  coarse-grain  poten¬ 
tial  parameters 
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3.  Fitting  of  atomistic  results  to  coarse-grain  potential  parameters  via  energy 
equivalence  and  consistent  mechanical  behavior 

4.  Validation  of  coarse-grain  model  with  full  atomistic  and/or  empirical 
results 

We  proceed  to  discuss  each  in  detail. 

2.3.1  Characterize  the  System:  Coarse-Grain 
Potential  Type  and  Quantity 

The  first  step  in  the  development  of  a  coarse-grain  model  is  to  determine  the  neces¬ 
sary  potentials  to  characterize  the  system.  In  general,  we  define  the  total  energy  of  a 
coarse-grain  system  as  the  sum  of  these  potentials: 


(2.12) 


E, 


'System 


where  <j)CG  is  the  defined  coarse-grain  potentials.  Although  a  simple  statement,  it  is 
not  trivial  in  implementation.  Indeed,  the  type  of  developed  potentials  can  serve  to 
both  broaden  and  restrict  the  applications  of  the  model.  Specifically,  the  coarse-grain 
potentials  must  be  complex  enough  to  represent  the  molecule  in  the  intended  meso¬ 
scopic  simulation.  One  must  consider  whether  provisions  are  necessary  for  fracture, 
intermolecular  interactions,  plasticity,  formation  of  secondary  structures,  or  a  host 
of  other  molecular  behaviors  depending  on  the  intent  of  the  investigation.  It  is  again 
stressed  that  a  coarse-grain  model  for  a  system  is  not  universal.  The  goal  is  to  utilize 
the  fewest  and  simplest  potentials  as  possible  that  represent  the  system  structure(s), 
mechanical  properties,  and  interactions.  Such  an  approach  can  be  facilitated  by  the 
utilization  of  harmonic  spring  potentials  to  reflect  mechanical  response,  where 


(2.13) 


Here,  refers  to  a  harmonic  spring  stiffness  while  i//  typically  refers  to  either  an 
interatomic  distance,  r,  or  an  angle,  6,  for  stretching  and  bending  deformations, 
respectively.  The  harmonic  spring  potential  results  in  a  linear  relationship  between 
force  and  deformation,  allowing  efficient  computation  of  interactions.  Indeed,  more 
complex  nonlinear  behavior  can  be  approximated  by  combinations  of  linear  func¬ 
tions  (such  as  bilinear  or  trilinear  formulations)  to  maintain  computational  efficiency 
in  lieu  of  the  introduction  of  more  complex  potentials  (such  an  approach  is  discussed 
further  in  the  ensuing  chapter).  A  possible  deficiency  of  harmonic  potentials  is  the 
continuous  increase  in  force  with  deformation.  Systems  subject  to  large  deforma¬ 
tions  will  subsequently  deviate  from  true  behavior,  and  thus,  without  provisions  such 
as  a  potential  cutoff,  the  use  of  the  harmonic  potential  should  be  limited  to  small 
deformation.  Selection  of  the  harmonic  spring  potential  can  also  be  affected  by  the 
atomistic  response  by  which  it  is  parameterized.  For  example,  a  linear  elastic  strain 
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response  can  be  sufficiently  modeled  by  harmonic  spring  potentials,  albeit  limited 
to  small  deformation. 

2.3.2  Full  Atomistic  Test  Suite 

The  finer-trains-coarser  approach  necessitates  the  parameterization  of  coarse-grain 
potentials  from  full  atomistic  results.  An  atomistic  test  suite  is  thus  developed  to 
obtain  the  necessary  mechanical  response  and  molecular  interactions  to  be  inte¬ 
grated  into  coarse-grain  potentials.  Typically,  a  relatively  simple  atomistic  simula¬ 
tion  is  applied  to  isolate  a  single  molecular  behavior,  mimicking  simple  material 
specimen  tests.  For  example,  uniaxial  stretching  can  be  applied  to  obtain  the  force- 
displacement  or  stress-strain  response  of  a  macromolecule,  allowing  the  calcula¬ 
tion  of  Young’s  modulus,  and  the  parameterization  of  coarse-grain  bond  strength. 
Further,  a  three-point  bending  test  can  be  utilized  to  determine  the  bending  stiffness 
of  a  molecule,  thereby  allowing  the  parameterization  of  the  coarse-grain  bending 
potential.  For  molecular  interactions,  pairs  of  molecules  are  simulated  to  deter¬ 
mine  relative  adhesion.  In  essence,  a  single  atomistic  investigation  is  implemented 
to  characterize  a  single  coarse-grain  potential,  ensuring  accurate  representation  of 
each  behavior.  The  number  of  required  tests  depends  on  the  number  of  coarse-grain 
potentials  (and  associated  parameters)  implemented  to  describe  the  system. 

2.3.3  Fitting  Coarse-Grain  Potentials 

The  fundamental  principle  underlying  the  parameterization  of  coarse-grain  poten¬ 
tials  is  the  assertion  of  energy  conservation.  Energy  equivalence  is  imposed  between 
potential  energy  functions  and  relevant  atomistic  energy  results.  The  basis  of  the 
formulation  of  the  energy  functions,  either  through  elastic  strain  energy,  deformation 
energy,  adhesion  energy,  or  other  techniques,  essentially  defines  the  accuracy  and 
behavior  of  the  representative  coarse-grain  system.  Two  such  approaches  of  energy 
conservation  warrant  separate  discussions:  (1)  the  application  of  direct  energy  equiv¬ 
alence  between  atomistic  and  coarse-grain  representations  and  (2)  the  assertion  of 
energy  conservation  via  consistent  mechanical  behavior. 

2.3.4  Direct  Energy  Equivalence 

Nonbonded  interactions  representing  macromolecular  adhesion,  attraction,  or  repul¬ 
sion  innately  have  no  mechanical  property  analogue  for  parameterization.  For  such 
energetic  potentials  as  the  Lennard-Jones  12:6  function,  or  Coulombic  interac¬ 
tions,  the  coarse-grain  equivalent  is  extracted  directly  from  the  energetic  results  of 
a  full  atomistic  representation.  Potential  energy  as  a  function  of  separation  can  be 
derived  between  individual  macromolecules.  The  coarse-grain  potential  can  then  be 
expressed  as 


(2.14) 
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where  (j'co  is  the  coarse-grain  potential,  R  is  the  macromolecular  separation,  and 
^  (Oifiij  ( r;/ )  is  a  weighted  summation  of  the  full  atomistic  potentials,  r/r  .  This  is 


essentially  how  coarse-grain  potentials  are  developed  for  a  small  number  of  atoms 
per  superatom  (such  as  MARTINI-type  force  fields),  where  the  superposition  of  mul¬ 
tibody  effects  is  apparent  in  the  mapping.  The  relation  becomes  more  complex  as  the 
number  of  incorporated/mapped  atoms  increase,  requiring  explicit  atomistic  simula¬ 
tion  of  the  desired  macromolecules.  The  above  relation  should  not  be  thought  of  as  a 
simple  summation  of  potentials,  or: 
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The  coarse-grain  potential  incorporates  effects  of  the  entire  represented  system 
(such  as  hydrophobic  and  other  solvent  effects,  electrostatic  interactions,  hydrogen 
bonding  between  molecules,  and/or  entropic  effects).  We  thereby  introduce  unknown 
weighting  coefficients,  (%,  in  the  summation  of  Equation  2.14.  The  complexity  neces¬ 
sitates  a  full  atomistic  investigation  to  derive  the  effective  coarse-grain  potential 
parameters. 

2. 3.4.1  Consistent  Mechanical  Behavior 

Bonded  and  angle  potentials  effectively  represent  molecular  mechanical  behavior. 
As  there  can  be  many  contributions  to  even  simple  mechanical  processes  (such  as  the 
breaking  of  hydrogen  bonds  or  solvent  friction  during  molecular  stretching  or  bend¬ 
ing),  we  introduce  conservation  of  energy  for  the  relevant  mechanical  response  (as 
opposed  to  the  potential  energy  of  direct  energy  equivalence)  resulting  in  consistent 
mechanical  behavior  between  full  atomistic  and  coarse-grain  representations.  The 
approach  is  to  apply  energy  equivalence  between  the  coarse-grain  potential  and  the 
observed  strain  or  deformation  energy  of  the  full  atomistic  system: 

(pCG(y/)  =  C%),  (2.16) 

where  (/>CG(y/)  is  again  the  coarse-grain  potential,  and  U(i//)  is  the  representative 
energy  function  for  the  mechanical  response.  It  is  an  underlying  assumption  of  the 
model  that  the  response  of  the  full  atomistic  system  is  accurately  described  by  the 
energy  function,  U(y/).  The  choice  of  the  strain  or  deformation  energy  function  is 
dependent  on  the  system  to  be  coarse-grained,  and  assumptions  of  elastic/plastic 
behavior,  material  isotropy,  system  failure  such  as  yielding  or  fracture,  etc.  affect 
the  definition  (and  interpretation)  of  atomistic  results.  As  an  example,  for  an  elastic- 
isotropic  material,  we  can  define  strain  energy  as 


U(e)  =  |J[<T][e]dV. 


(2.17) 


30 


Multiscale  Modeling:  From  Atoms  to  Devices 


To  equate  with  we  must  formulate  U  as  a  function  of  i//.  For  axial  stretching, 

this  requires  the  formulation  of  strain  as  a  function  of  bond  length,  e(r),  while  for 
shearing  this  requires  the  formulation  of  strain  as  a  function  of  shearing  angle,  e(9). 

Individual  parameterization  of  coarse-grain  potentials  via  consistent  mechanical 
behavior  has  an  inherent  limitation  by  the  decoupling  of  mechanical  behavior.  To 
illustrate,  as  bond  and  angle  coarse-grain  potentials  are  independent  under  simula¬ 
tion,  there  is  no  coupling  between  stretching  and  bending  of  the  molecule  and  effects 
such  as  twisting  or  torque  are  completely  neglected.  However,  under  full  atomis¬ 
tic  representations,  such  coupling  can  serve  to  increase  stresses  and  induce  failure. 
Ultimately,  the  model  is  tailored  to  the  desired  application,  and  systems  with  obvious 
or  significant  coupling  effects  should  be  avoided. 

2.3.5  Validation 

Once  the  coarse-grain  potentials  are  developed,  validation  of  the  coarse-grain  model 
is  a  necessary  step  to  assure  an  accurate  representation.  The  primary  basis  for  vali¬ 
dation  is  through  a  direct  comparison  to  full  atomistic  simulations.  Such  an  approach 
may  seem  counterintuitive,  redundant,  or  self-serving  as  the  coarse-grain  potentials 
were  developed  from  such  simulation  results.  However,  the  test  suites  applied  were 
intended  to  result  in  a  single  behavior/response  for  each  coarse-grain  potential.  The 
full  coarse-grain  representation  can  possibly  involve  the  interaction  and/or  coupling 
of  individual  potentials.  The  combination  of  two  behaviors  (bending  and  stretching, 
for  example)  must  be  explored  to  justify  the  coarse-grain  model.  Relatively  simple 
validation  simulations  can  be  designed,  combining  multiple  system  behaviors,  and  a 
one-to-one  comparison  made  between  full  atomistic  and  coarse-grain  results. 

Due  to  limitations  of  full  atomistic  simulations,  such  validation  is  restricted  to 
component-level  behavior  in  which  the  coarse-grain  model  is  specifically  developed 
to  circumvent.  If  available,  secondary  support  and  model  validation  is  found  in  exper¬ 
imental  data.  Experimental  techniques  such  as  nanoindentation  or  atomic  or  chemi¬ 
cal  force  microscopy  (AFM/CFM)  can  directly  probe  materials  at  the  mesoscale  and 
determine  system  level  characteristics  and  mechanical  properties.  The  development 
of  coarse-grain  models  is  intended  to  elucidate  such  system-level  behavior  and  thus 
correlation  with  experimental  results  is  essential.  Reciprocally,  an  accurate  coarse- 
grain  model  can  serve  to  validate  and  support  experimental  results. 

It  is  noted  that  advances  in  experimental  techniques  are  continually  probing  smaller 
scales,  reaching  atomistic  precision,  and  allowing  single  molecule  investigations  [such 
as  optical  tweezer  methods,  atomic  force  microscopy  (AFM),  etc.].  Such  data  can  serve 
to  reinforce  both  the  full  atomistic  representation  and  the  developed  coarse-grain  model, 
thereby  increasing  confidence  in  both  representations  and  supporting  the  finer-trains- 
coarser  approach.  Figure  2.6  provides  a  schematic  of  the  relationship  between  atomistic 
and  coarse-grain  representations  with  experimental  data  for  model  validation. 

2.4  SUMMARY  AND  CONCLUSIONS 

The  interplay  between  material  properties  and  structural  arrangement  of  hierar¬ 
chical  systems  provides  unique  and  robust  mechanical  behavior  that  transcends 
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FIGURE  2.6  Validation  of  coarse-grain  model  through  component  behavior  at  atomistic 
scale  and  system  characterization  from  experimental  data.  Experimental  data  ranges  from 
component  (molecular)  to  system  (mesoscopic)  length  scales,  providing  a  reciprocal  loop  of 
validation  with  both  atomistic  and  coarse-grain  models. 


multiple  length  scales.  Coarse-grain  models  provide  an  efficient  means  to  simulate 
and  investigate  such  systems  in  which  the  desired  behavior,  property,  or  response  is 
inherently  at  the  mesoscale.  The  relatively  few  parameters  and  potentials  that  define 
a  coarse-grain  representation  must  accurately  represent  the  full  atomistic  behav¬ 
ior  of  the  structure  in  both  mechanical  response  and  intermolecular  interactions. 
Concurrently,  the  diverse  nature  and  possible  components  of  hierarchical  systems 
makes  a  general  formulation  difficult.  We  thus  introduce  and  discuss  a  system- 
dependent  finer-trains-coarser  approach  to  coarse-grain  model  development,  requir¬ 
ing  the  parameterization  of  coarse-grain  potentials  via  full  atomistic  results  in  a 
consistent  and  systematic  manner. 

A  well-defined  coarse-grain  representation  simultaneously  allows  investigation 
of  molecular  structures  and  behavior  at  multiple  levels  of  hierarchy  while  intrinsi¬ 
cally  maintaining  full  atomistic  behavior.  This  simplification  also  allows  a  pragmatic 
approach  to  system  analysis,  delineating  global  system  behavior  from  constituent 
element  properties,  thereby  encompassing  multiscale  structure-property  relations. 
Extension  of  accessible  time  and  length  scales  can  allow  direct  correlations  with 
empirical  investigations,  providing  a  novel  tool  for  experimental  design  and  material 
characterization. 

Coarse-grain  model  formulation  can  be  described  in  terms  of  four  guiding  prin¬ 
ciples:  (1)  selection  and  characterization  of  governing  potential  types  and  quantity; 
(2)  application  of  full  atomistic  test  suite  to  reflect  independent  behaviors  and  deter¬ 
mine  relevant  coarse-grain  potential  parameters;  (3)  fitting  of  atomistic  results  to 
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necessary  potential  parameters  via  direct  energy  equivalence  of  conservation  of 
energy  through  consistent  mechanical  behavior;  and  (4)  validation  of  the  developed 
model  via  comparison  with  full  atomistic  results  (component-level)  or  correlation 
with  empirical  data  (system-level). 
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